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Introduction 
In this paper we develop the theory of equations of the form 
(0.1) a{x)i)(T{x)) + (3(x)ip(x) + -/(x)^(t^ 1 (x)) = \ip(x) 

generated by a bijection r : X — > X of the real line subset X C R. 

The equations of this type can be regarded as an alternate discretization of the 
second order differential equations including Schrodinger equation and generaliza- 
tion of difference and g-difference equations. They also emerge from the change of 
variables in difference equations. On the other hand the functional equations by 
themselves are of interest and have many important applications. They have been 
investigated in many monographs and papers from the point of view of functional 
analysis, especially C*-algebras methods, see e.g. [21E]. We hope that our approach 
can be applied also in the numerical analysis of some problems. 

Our approach is based on r-difference and r-integral calculus, see e.g. ^T], which 
is a direct generalization of the standard one and is ideologically related to the 
analysis on time scales. This gives us a possibility to generalize the factorization 
method elaborated effectively by many authors in the differential case, e.g. see 
El- Our approach allows one to consider this method in q-case, see jH] for 
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application, and due to generality of r has the benefit that it is equivariant under 
the change of variables. 

Section H defines basic notions and operators of r-calculus as well as formulates 
some relations between them. In Sectional we introduce and study Hilbert spaces 
with the scalar product given by r-integral on the orbits of the semigroup of nat- 
ural numbers (N U {0}, +) or group of integer numbers (Z, +) acting on X by the 
iterations of the bijection r. Then, in Section we construct the chain of eigen- 
problems in these Hilbert spaces by the application of the factorization method. 
This construction gives a possibility to solve the eigenproblems for the chain of 
Hamiltonians if one knows a solution to the initial one. 

In Section we discuss the change of variables and induced transformation of 
Hilbert spaces, eigenproblems and parametrizing functions. Some remarks about 
equivalence of one r-difference operator to another are also presented. 

In Section|3we show that our method includes the one investigated in the theory 
of g-Hahn orthogonal polynomials, which are q-deformation of the classical ones, 
see ^2]. In this section we also present some examples. 

Some general facts about r-equations is presented in Appendix A, i.e. some for- 
mulas for solutions. Reader can find there easy results which are the generalization 
of theorems about g-Riccati equation from • 

Let us stress finally the importance of the functional equation 10. Ill in the theory 
of discrete and g-discrete polynomials, see ^2]. This link is also important to the 
integration of some quantum optical systems, see [2]. 

1. Calculus generated by the bijection of subset of the real numbers 

Let us introduce preliminary definitions and facts of the calculus generated by 
the bijection t : X — » X of the subset Id (see also |1L). Let 3-{X) denotes the 
algebra of all complex valued functions ip : X — ► C. Let O(x) = {T Il (a;)|n G Z} be 
the orbit of point x under action of group Z given by r and + {x) = {r n (x)\n <E 
N U {0}} — orbit under action of semigroup N U {0}. 

The subject of the paper is linked to the investigation of the following linear 
operators 

(1.1) Tip(x) := <p(r(x)) 



(1.2) Mfip{x) := f(x)<p(x) 

ip(x) - (p{r(x)) 



(1.3) d T ip{x) 



x — t{x) 



where / £ T{X). 

The operator d T is a direct generalization of g-derivative d q , which corresponds 
to the special case when t(x) = qx. We shall call it r-derivative. It satisfies the 
following analogue of Leibniz rule 

(1.4) d T {ipip) = (T(p)d T tp + i;d T (p. 
We define r-integral as the solution of the equation 

(1.5) d T ip = V> 
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for given ip S F{X). Applying the operator T n to both sides of lll.rijl one obtains 

(1.6) V{r n {x)) ~ <p(r n+1 (x)) = (r n (x) - r n+l (x))iP(r n (x)) 
for n = 0, 1, 2, . . ., where r n = tq. . .qt and thus 

n times 

oo 

(1.7) p(x) - ^(r°°(x)) = £> n (x) - r^ 1 ^))^^^)). 

ra=0 

The formula l|1.7|) has sense only if t"(x) has the limit when n — > oo and (/? is 
continuous at r°°(a;) £ X. Hence, if these conditions are satisfied, one can define 
the r-integral: 

(1.8) / ij(t)drt = / mdrt := Y,(T n (x) - T n+1 ( X ))^(T n (x)). 
Jr°°{x) JO+{x) n=Q 

Let [a, b] C I be an interval such that r oc '(a) = t°°(6) then we shall define the 
integral of ip over [a, b] by: 

pb p pb pa 

(1.9) / ip(t)d T t = I ip(t)d T t := / V(*)<M - / ip{t)d T t. 

Ja J 0+ (a)UO+ (6) Jr=°(b) Jr°°(a) 

Let us note that integral is only formally over [a, 6] and it is based on + (a)LiO + (b). 
Similarly we will define r-integral over 0(x) by 

p oo 

(1.10) / mdrt:= (r n (x)-r n+l (x))rP(r n (x)). 

The terminology introduced above is justified by the following properties of the 
r-integral: 

(1.11) [ (d T ip)(t)d T t = i;(b)-i>(a), 



(1.12) d T / 1p(t)d T t = 1p(x). 

Jt°°{x) 

Moreover we have the following formula for the change of variables 

(1.13) / (Tip)(t)p(t)d T t = / ^idrT^miT- 1 p){t)d T t. 

J a J r(a) 

Let us denote by T>{X) C F{X) the subalgebra of functions with finite support. 
This subalgebra is preserved by T and d T , as well as by Mf. 

For affine bijection r(x) = qx + h all formulae above reduce to the corresponding 
formulae of g-difference (h = 0) or /i-difference (q = 1) calculus. If q —> 1 (h = 0) 
or /i — > (g = 1) one obtains the formulae of the standard differential and integral 
calculus. Formulae (11.1 If) and (11.1 2|i become then the fundamental theorem of 
calculus and l|1.13fl reduces to identity. 
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2. HlLBERT SPACES WITH SCALAR PRODUCT GIVEN BY T- INTEGRAL 

Let us fix some weight function p : X — > R. We can define the scalar product by 
the r-integral over the set S by the formula 



where S is one of the following sets O(xo), + (a), + (a) U + (b) and ip,ip are 
complex valued functions defined on X. 

For the positivity of <EjJ we shall assume that (r™(6) - r n+1 (b)) p(r n (b)) ^ 
and (T n (a) — r n+1 (a)) p(r n (a)) ^ 0, n £ N U {0} for semigroup orbits case or 
(T n (xo) — r n+1 (xq)) p(r n (xq)) 5* for group orbit case. 

Let us denote by L 2 (X,pd T ) the Hilbert space of square-integrable functions in 
the sense of <2.U . It is clear that T>(X) C L 2 (X,pd T ) for any weight function p. 
One shall keep in mind that elements of this space are classes of equivalence of 
functions from !F{X), which can be identified with functions on S. 

If we consider semigroup orbit case + (a) U + (b) and let b — t n (a) for some 
N e N, then L 2 (X,pd T ) is finite dimensional. This case will not be discussed in 
the paper. So, it will be assumed everywhere that b ^ t n (a) for all N E N. 

If one admits the case p(x) — in 112. If) for some x € S, then the Hilbert space 
reduces to direct sum of Hilbert spaces realized by the functions on orbits or finite 
sets. So, one comes back to the cases defined above. The case of two semigroup 
orbits is also direct sum and it is considered here only due to the analogy with 
g-Hahn polynomials. 

Let us now consider the operator T acting in Hilbert space L 2 (X,pd T ) defined 
on V{X) by fTTTjt . Note that V{X) is dense in L 2 (X,pd T ). 

Proposition 2.1. The adjoint operator T* is the weighted shift operator given by: 



The domain T>t* of T* contains V(X). 

Let us mention the following properties of the operators T and T* . For the group 
orbit case one has identity 



(2.1) 




(2.2a) 





(2.4a) T*Tip = pip. 

For the semigroup orbits case one obtains 

(2.4b) r rT<p = f i(l- X {a,b})<P, 

where X{a,b} ls the characteristic function of the subset {a, b} C X and thus M X{a b) 
is a projector onto this subset. For both group and semigroup orbits cases one has 

(2.4c) TT*ip= (poT)ip. 
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The operator T is bounded if and only if sup \p(x)\ < oo and then ||T|| = ||T*|| = 

xes 

(sup\n(x)\)i. 
xes 

Let \a and \b be the characteristic functions (indicators) of + (a) and + {b). 
The orthogonal projectors M Xa and M Xb map T>(X) into itself. Since operator T 
preserves subspaces L 2 (0 + (a) , pd T ) and L 2 (0 + (b) , pd T ) , it can be reduced to the 
components of the partition 

(2.5) L 2 {X lP d T ) = L 2 (0+(a),pd T )®L 2 (0 + {b),pd T ) 

in the sense of pQ. Thus we can consider each orbit separately. 

Let us now consider the sequence of Hilbert spaces Tik '■= L 2 (X, Pkd T ) with the 
weight functions pk, k E N U {0} related by the recurrences 

(2.6) pk+i = VkPk 
and 

(2.7) p k+ i = T{B kP k) 

for Bk,i]k £ Since we demand two distinct formulae for pk+i, we need to 

add the consistency condition on T)k and Bk, namely 

(2.8) T(B k pk) = VkPk 

and in semigroup case we additionally impose a boundary conditions 

(2.9) Bk(a)p k (a) = B k (b)p k (b) = 

for any k £ NU {0}. 

Let us denote by gk the ratio of Bk and Bk+i- 

(2.10) B k +i=g k Bk. 

Then we can easily see that formulas l|2.(ill - lj2.H|l are automatically satisfied if (12.81 — 
(I2.9JI are valid for k = and 

(2.11) rjk+i = T(gkVk) 

for any k £ NU{0}. 

By introducing the function 

. / n Bk(x) — rihix) 
2.12 A k {x) = kK ' lkK } 

X — T(X) 

and expressing r/k by Ak and Bk, one can rewrite equation l|2.8Jl in the form 

(2.13) d T (BkPk) = A k pk 

In the case of r(a;) = qx and when q — ► 1, the equation (12 . 13f) corresponds to 
the well known Pearson equation of the theory of classical orthogonal polynomials. 
This motivates us to call the consistency equation 112, 8J1 the r-Pearson equation. 

We shall now consider the operators defined formally by ljl.ljl - ljl.3|l as operators 
in Hilbert spaces. For / e J-(X) we shall consider the operators of multiplication 
by the function / as 

(2.14) Mf.Hk^Hk+i. 
The operator T acts in each of the spaces separately 

(2.15) T-.Uk^Uk 
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and the r-derivative becomes now the operator 

(2.16) d T = M (id _ T) -i(l -T):H k ^ H k +i- 

Since operators M/, T and d T are in general unbounded, we take V(X) as their 
domain and then close them. We will denote their closure by the same symbols. We 
also omit indices k in symbols of operators M, T and d T to simplify the notation. 

In these settings we can apply equation 112. 8J1 to calculate explicitly the term \ik 
appearing in formula 112.311 for T*: 

(2-17) f i k (x)=d T (r- 1 )(x) fff„ ■ 

Proposition 2.2. 

The adjoint operators MJ : Hk+i — > T~tk and d* : Tik+i — > 7~tk, are given by 

(2.18) m; = M ml 

(2.19) d; = {l-T*)M m{id _ T) ^. 
The domains of Mj and d* contain T>(X). 

Proof: By means of i'l.CA) and <2.H we get 



(/¥#)k+i = / f(x)ip(x)il>(x)pk+i(x)d T x 



s 



= / (p(x)ip(x)f(x)r) k (x)pk(x)d T x = {<p\fr)ki>)k 



s 

for 99,^/) € "D(S), Formula for 9* follows directly from Q2.18I1 and definition of d T 

$nfy . □ 

3. Factorization method 

The proposed method is the direct generalization of method for the second order 
differential equations, see 0E3|- Instead of the equation of the form 

(3.1) a{x)^{T{x)) + P{x)tp(x) + 7 (x)V>(t- 1 (x)) = Xtp(x), 

where A £ C and a, f3, 7 £ J~(X), we will consider the sequence (chain) of the 
factorized equations 

(3.2) (A* k A k ^ k )(x) = XMx) 

on ijjk £ Hk with the operators A& : Hfc — > Hfc+i and their adjoints A£ : Hfc+i — > 
Wfc given by 

(3.3a) A k :=M hk d T + M fk 



(3.3b) A£ = (l-T*)M Wid _ r) _i +M w/feJ 

where : X — > K, /i^ : X — > R and fe £ N U {0} . We assume that is a closure of 
the operator defined by l|3.3all on T>(X) and A£ is automatically closed. Then by 
von Neumann theorem (see [Q) A£A& is self-adjoint and thus closed. The essence 
of the approach is to postulate the relation 

(3.4) A k A* k = d k Al +1 A k+ i + c fe , 
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where c/-, dk S C and t £ N U {0}. It follows from I|3.4J1 that given a solution 
^fc G Wfe of equation i|3.2|) , the function 

(3.5) Vfc+i := A fe ^ fc e W fe+1 

is a solution of (13 .21 for fc + 1 with A^+i = 7° k ■ So, starting from the equation 
fl.-S.2ll for k — and its solution f/'o, one generates the family of equations fl.-l.2ll with 
solutions given by 

fe-i 

(3.6) i> k = [] AiVo- 

The operators and A£ depend on the functions f k , /ifc as well as on n k and 
Bfc related by the recurrences 12. 1H and 112.1011 respectively. 
Let us introduce the function 

(3.7) ^ k{x) . = h{x) + J}±W 

X — T(X) 

which we will use instead of f k . 

Proposition 3.1. Condition 113.411 is equivalent to the transformation rule 

(3.8) ipk+ihk+i = 3^ T[ — 

d k V 9k 

and the system of equations 



d k g k (x)B k (x)(h k +i(T 1 (x))) 
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(<p k (x)) Vk(x) + c k 



(x — t(x))(t 1 (x) — x) 
ro^ 1 ( d k g k {T{x))B k {T{x)){h k {x)) 2 2 {h k {x)f 

tS^R) { (r(x)-r»(x))(x-r(x)) " ^^(W^)?. 
on i/ie functions B k , r\ k , ip k , g k , h k , hk+i and constants c k , d k , where k £ NU{0}. 

Proof: We compute explicitly left and right hand sides of l.'t.4t . By substituting 
the IJ3.3I) . to (I3.4|l and applying l|2.4hjl one has 

A fc A* = -^(^kiTjx^nkiTjx)) T {h k {x)) 2 B k {r{x)) 



x — t(x) ( t ( x ) — T 2 (x))(x — t(x)) 

(3.10 + ^ fc x ) r/fc x rl T 

x — t(x) 

and 

a * a = -yfc+i(a:)?7fc+i(a:)fefc + i(x) j9fc + i(x)(/i fc+ i(r~ 1 (x))) 2 
fe+1 rc-rfc) (a;-T(a;))(T- 1 (a;) -2) 

(3.11 +(Vfc+i(a:)) %+i(ai) -r-. T . 

X — T(X) 



In (13 . X If) we have removed the projector 1 — X{a.b} appearing in (I2.4bll for the case 
of orbits of semigroups. It was multiplied by a function B k+ i, which by Ij2.9|l is 
equal to zero at the points a and b anyway. 

Substituting (|3 . lOfl and l|3.1H into (13.411 and comparing coefficients in front of 
the operators T, id and T" 1 , one obtains 

(3.12) hk(x)(p k (r(x)r)k(T(x)) = d k hk+i(x)ipk+i(x)r) k +i(x), 
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Bfc(r( rEy ) i^ +(^(^)) 2 %w 



(t(x) — r 2 (a;))(x — t(x)) 

(313) =4 ( (,-r(,))(r-H,)-,) 
and 

(3.14) y h {x)B{x)h k {T- l {x)) = dfcSfc+i^^fe+^r-^a;))^!^- 1 ^)). 

Using the transformation rules 112.1111 and 12. Kit one sees that conditions H3.12JI 
and Q3.14I1 are equivalent. They give the transformation rule l|3.8|l 



(3.15) ipk+ihk+i = h k T ipk 

\dkgk 

and condition l|3.13|l can be written in the form l|3.9|) . Domains of AfcA£ and 
A^ +1 Afe + i coincide due to closeness. □ 

System of equations II3.9J1 can be regarded as an equation for gauge functions gk 
(defined by H2.11l) - H2.1f)J) l or transformation rule for hk given in non explicit way. 
Check Example 15.11 for explicit formulae for gk in the case hk = 1 and Cfc = as 
well as convergence conditions. In the case t(x) — qx and gk = const solutions to 
113. 9JI have been classified in 

If we have the solution of l|3.9ll (both hk+i and gk) one can express the func- 
tions (Bk+i,rik+i, fk+i) of the (k + l) s * step in the terms of the previous functions 
(B k ,VkJk)- 

Functions Bq, 770 , <^o and Ao are related to initial equation l|3.1H by 
(3.16 a{x) = — , 

X — T\X) 

(3.17) 0(x) = (Mx)) 2 Vo(x) + 



{x — t{x))(t 1 {x) — x) 

(3.18 tW = n , 

X — T(X) 

and 

(3.19) A = A . 

The formulae inverse to (|3.1fi|) - (l3.19|) have the form 

„ , s (x — t(x))(t (x) — x) ( ' „ , , ,<j9n(x) 

3.20 B (x) = ^ \. \ > 1 (3(x) + (x - r (x a (x )^j-L 

(x — r(x))a(x) 



(3.21) r) (x) 



ipo{x)h (x) 



where ipo and ho have to be chosen in such way that their ratio satisfies the equation 
of the r-Riccati form l|A.llj) : 



^ o(r( - T)) ^a{r{x))(r{x)-T^x)) 
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By this approach we have reduced the problem of construction of the chain of 
Hamiltonians to finding the solution of Ij3.9l) and choice of initial conditions is ruled 
by 

4. Change of variables 
We are going to consider an invertible change of variables 

(4.1) k:X — ► Y. 

It is convenient to assume that k is homeomorphism in order to preserve the limits. 
Such transform allows us to modify the function r in considered equations. Let 
us note that in differential equation change of variables keeps the differentiation 
operators and only changes the direction of derivation and multiplies them by some 
function factor. In discrete case one obtains in fact different derivations. 
Change of variables defines the map from T{Y) to F{X) by 

(4.2) K :F(Y)3f^Kf = fo K er(X). 
We can associate to T the operator 

(4.3) f = K~ 1 TK 
acting as the shift by 

(4.4) T=KOTO K~ 1 

on J-(Y). One calls functions r and f equivalent if they are related by \AAl for 
some K. Operator Mf, f € F{X), transforms to 

(4.5) K~ l MfK = Mf 0K -i. 

Thus similarity transformation K~ x ■ K maps shift operators to shift operators and 
multiplication operators to multiplication operators. However it turns out that 
derivation is preserved by this transform only up to the function factor: 

(4.6) K-'drK = M (dTK)oK -id f . 

For example let X = [0, oo), t{x) — qx, Y = R for < q 1. Then k(x) = ln(x) 
transforms r to 

(4.7) f(y)=y + \nq 
and 

(4.8) KdnR- 1 = M -m q d f . 

Thus we can map g-difference equation to /i-difference equation although there 
may appear function factor changing the form of equation. 

For given (X, r) and (Y, f) in general it is not possible to construct such k that 
(I4.4|l would be valid. It is the case for example for X = Y = [0, 1], t(x) = x 2 and 
f (y) — |y 3 — |y 2 + |y. It follows from the fact that k maps fixed points of r into 
fixed points of f. Since t has two fixed points (0 and 1) while f has three (0, \ 
and 1), thus f is not equivalent to r. In some cases however one can restrict X and 
Y in such way that the function would become equivalent but it may happen that 
K is not given in terms of elementary functions thus rending further calculations 
extremely complicated. 

Let us now consider change of variables defined by H4.1H from the point of view 
of Hilbert spaces and introduced operators and functions. Let us consider Hilbert 
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spaces Hk of functions on the set Y denned in such way that it makes K an unitary 
operator mapping Hk to Hk- Simple calculations shows that we have to consider 
the spaces H — L 2 (k(X), p k df) with weight given by 

(4.9) pk^idfK- 1 ) R- l Pk . 

In order to keep formulas for transformation of pk the same as \2M and H2.7J1 one 
has to define 

(4.10) B k = -^-l K ^B k 
and 

(4.11) i% = if~V 

These functions automatically satisfy f-Pearson equation and obey the same trans- 
formation laws with function 

(4.12) g k = K-'gk. 

If we put the following transformation rules for functions fk and hk 

(4.13) h k = K- l (h k d T {K)) 



(4.14) f k = K-'fk, 

then operators A& = K A k K and A£ = K~ 1 A* k K are of the same form as 
original with r replaced by f. Clearly A* k is an adjoint of Ak since K is an unitary. 
Let us note that due to l|4.13ll the condition h k = 1 is not preserved by the change 
of variables except in the case of translation k(x) = x + c, c e K. 

Solutions of considered second order equation 13.20 transform as follows 

(4.15) i> k = A-ty*. 

Thus we were able to map the chain of Hamiltonians on X to the chain of 
Hamiltonians on Y preserving all relations introduced in Section and 03 This 
shows that considered method is equivariant with respect to the change of variables. 
It was not true for the case of g-difference equations. 

5. Some examples 

We are going to present several examples for the application of factorization 
method described above. 

5.1. Particular solution of 13.9J1 . We are going to construct a solution equation 
13.9|l in terms of r-integrals by application of propositions from Appendix lAl 
We consider the simplest case hk = 1- Let us introduce the function 

, K1 s e , s f / \ \ 2 / \ d k g k {x)B k {x) 

(5.1) := Vk(x) - {x _ T{x)){T - 1{x) _ x y 

After substituting £k instead of gk the equation 13.90 takes the form 

( (rM-^tS-rW) -Ck)Ur(x))+c k (Mr(x))) 2 Vk(T(x)) 
(5 ' 2) (<fik(r(xWVk(r(x)) - Z k (T(x)) 

which is also a r-Riccati equation lA.llll . It is more convenient than the form 13.9|l 
because in the case of Ck = one can apply the Proposition IA.1I from Appendix 
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El which allows one to write down formulae for solutions. To that end we have to 
invert the fractional map on right-hand side of (15. 2J1 and use it to define matrix 
Afe(x) appearing in l|A.9j) . We get: 

A (x-t(x))(t(x)-t 2 (x)) \ 

(5.3) Afe(x) = (x-T(x))(r(x)-J(x)) Vk (r(x)) 2 ri k (r(x)) ' 

\ U B k (r(x)) ) 

Since this matrix can be singular at t°°(x), we perform r-Darboux transform 1|A,34J| 
to Ak(x) with 62 — A and Si = 0. In order to apply Proposition IA.3I we have to 
compute A' fc defined by IIA.8J) : 

/q (t(x)-t\x)){x-t™{x))>- \ 

(5-4) A k (x) = I y,^^)) 2 ^^^ (t^-t^Xx-t""^)^ • 

\ U x-t(s) B*(t(x)) (t(x)-t<*>(x))A / 

Let us put A = — 1. If the functions ft, r\k and B k are regular and non- vanishing 
at t°°(x) and r is differentiable at t°°(x) then we have 

(5.5) lim A^«(*)) = f° - ~ (^M))^ . 

n-oo / 

In such a way the singularity in A fc is removed. If we assume that r satisfies the 
condition i) of Proposition IA.3I then the resolvent matrix is well defined and by 
Proposition I A . II we find 

(5.6) 

(I F k (x) 

KM = 

where Fk(x) is given by 



in 



(T(t)— r°°(t))(t— r(t))(T(t)— H (t))y fc (T(t))"„ fc (T(t)) 



ex P (r T i (x) - ( '-T-^ (T(,)) <M 



, (r(t)-r°°(t))(t-r(t))(r(t)-r^t))y fc (r(t))^ fc (r(t)) 

(.) = ex P 1 / (t -;°° (t) ;y- (t)) <m 1 x 

(r(t)-^(t)) x 



r oo (a!) (*-T~(t))B fc (r(t)) 



« In- 



/ r oo (x) s-r(s) 
Thus, by l|A.15jl and l|A.36jl one gets 

&(x) = (x-r~(x))- 1 exp -/ ^ ] x 

X (&£) 7 



t-l_f a: (T(t)— r 2 (t)) / ft ln [ r (,)- T - M )(,- T °°(,))(,-,;,))( TM - T H,)) ftM ,))a, t ( TM ) , \ , . 

SfeO J T =°(i) (t— r°°(i))B fc (r(t)) ^Ut"^) s-t(s) u T h \U t l 

From lO by application of iJEH, E3, (ZE3 , @IIH and we can 

write down explicit form of operators Afe + i and A£ +1 . We will not do it here due 
to length of formulas. 



B fc (r(s)) 
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5.2. g-Hahn polynomials. Let us show now the well-know application of the 
method described in a previous section to g-deformed classical polynomials. We 
consider t(x) — qx for < q < 1 and to stress it we replace the symbol d T and 
T by d q and Q respectively. We assume that Bo and A$ defined by 112.1211 are the 
second and first order polynomials respectively. We consider a function po satisfying 
g-Pearson equation H2.13I) and we put fo = 0, ho = 1. In this case the equation 
113. If) is known as Hahn equation. Its solutions are given by basic hypergeometric 
series, see |ESEj, and we know from the theory of ^-classical orthogonal polynomials 
that for certain choices of Ao , we obtain polynomials which belong to Ho and form 
orthogonal system with respect to weight function po, see |12j . 

We want to keep the same type of equation for all k, i.e. we require that all 
f k = 0, h k = 1. We shall find the conditions under which this requirement can be 
satisfied along with relation (13 .4|) . 

Equation 113. 8JI reduces itself to 

(5.9) gkdk^q^ 1 - 

The transformation rule for % 12. lit can be rewritten in terms of A k as follows 

(5.10) A k+1 (x) = qQA k (x)+d q B k+1 {x). 

Thus A k and B k remain a polynomials of degree one or two respectively for all k. 
We can express function defined by llo.lll in terms of our variables as 

(5.H) a- = tt^V 

(1 - q)x 

and equation ll."i.2ll is satisfied for c k — —d q A k , which is a constant for all k. Let 
us mention that since is singular we can obtain this solution from H5.2J1 by 
the product l|A.15j) described in Appendix El only by applying singular Darboux 
transform l|A.34j) with Si — 62 = 1- 

Formula ^ yields that if {P*} is an OPS for p k then {d q P%} is an OPS for 
Pk+i- We can reverse this procedure. Namely, by applying A£ to both sides of 
i|3,5Jl we gather that 

(5.12) Vfe = ^-AJVfc+i. 

This allows us to express any polynomial by the formula 

(5-13) P* = p k n A* k+1 ...Al +n l, 

where p^ ^ are some normalization constants, see |12j . 

We can generalize this result by dropping the requirement that f k = for k ^ 1 
and then put c k = 0. Then we can solve the equation H5.2H by the formula H5.8H . We 
can express that solution by basic hypergeometric function. This case is presented 
as an example in [J]. 

By the application of the change of variables described in Section 01 one can 
obtain solutions to another family of equations with t(x) ^ qx. For example, let 
k(x) = e x . Then f(y) = y q . Hahn equation transforms into 

(5.14) A£A fc ^ fe = 0, 

where 

("i ^ = ter^ 1 * + (^) 2 i a 4) + 



1 - q J Iny 

Its solutions are related to orthogonal polynomials by formula Q4.15I1 and are or- 
thogonal with respect to the scalar product with weight function 

(5.16) p k (y) = -^(1 - g)A^Vfc(y)- 

v - y q 

5.3. The case g k = const and t(x) = qx. This example is based on the paper 
We assume that 

(5.17) t{x) = qx for < q < 1, 

(5.18) h k = 1, 
(5-19) 9k=q~\ 

(5.20) d k = 1 

and that B Q is a constant function. From 112.101) we gather that 

(5.21) B k = q- 2k B . 

Let a k := (ifik) 2 Vk- Equation (13 . Of) now takes the form 

(5.22) a k (x) = q~ 2 a k {qx) - c k . 
By substituting into equation l|5.22l) the series 

(5.23) a k (x) =.T A ^v n , 

nGZ 

where A 6 [0, 1), we obtain 

(5.24) a k (x) = b k x- 2 - 

1 — <jr 

for any 6 C. In order to preserve the transformation rules H2.11I) and 113.811 we 
have to put b k — q 4k bo and c k — q 2k co. Under these assumptions the operators A k 
and A£ take the form 

(5.25) A fc = - 1 Q + q 2k (Q k vo), 

(i - q)x 



q 2k b x^ -j^ B, 



(5 ' 26) (QVo) q 2k (l-q)x^ ■ 

Thus we have solved the factorization problem for the operators 

AkAk { X Vk {x){l-qf{l + q) X <p k (x)(l-q))* 

3 (1-3) V 3(1-3) / (1-3 ) 

The relation 

(5.28) K-i^k = 
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implies that A* k A k tp k = (Afc_iA£_ a - Ck-i)ipk = -Ck-itpk- Thus ^ is an eigen- 
vector of AfeA^ with the eigenvalue —Ck-\. We easily find the solution to i!5.28|l 
by the use of IIA.3I) : 



In 



q)qta k {t) 



(5.29) ip k = nexpj — v t^l-q^ ~ dqt ' 

where K is a constant. Let us assume that H5.29J1 is convergent. We have to check 
if ipk G Wfe, i.e. if it is square integrable with weight function p k emerging from 
g-Pearson equation l|2.8|l . To that end we shall calculate function ip k p k . Let us 
rewrite Pearson equation in the form 

„ n x Q-k QPk 

(5.30) — T — = . 

From II5.28|I we conclude that 

0}h B k a k 

(5-31 —— = J- r-3 

Vfe (1 - q)q A x ip k 

and thus 

(5.32) rtpk = (b k Y^x 2 ) (1 ^ Q(V^)- 

We can rewrite this equation in the following form 



(5-33) nPk=l[l--p) QWP*). 



where /3 = y bfc ^ fc <? - and 7 = — — bfc . If we request that 7 = 1 by appropriate 
choice of b k and then we have the following solution 

(5.34) (M*))V(z) = m(-|;<?) (J;<? 

where 

(5.35) (a; q) n := (1 - a)(l - «a) • • • (1 - g n - x a) 

and /i € M\{0} is such that the scalar product is positively definite. Let us estimate 
the g-integral of this function. 

(5.36) J^{Mx)) 2 Pk(x)d q x = p{l-q)a^q n ^-^-q^j (^1?) 

Since (±a/(3;q) n is convergent when n^oofor0<(7<l then (±aq n //3;q) 00 
is necessarily convergent to 1. Thus it is bounded from below and above. Let us 
denote these lower bounds by m± and upper bounds by M±. Thus 

r-a 00 

/ {ipk(x)) 2 p k (x)d q x sc |^o|iriax(|M_| Jm-|)max(|M+|,|m+|)(l - g) ^ = 

J n=0 



(5.37) = |/za|max(|M_ , |m_|) max(|M+| , \m+\). 

So this integral is finite as well as (ip k \ip k ) k = f i^ k pkd q x and thus ip k £ H k 
any choice of limits of integration fulfilling 112. 9J1 and constants such that 7=1. 
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Thus by 113. 5 jl we have 

(5.38) = A k+N-i ■ ■■■■ Afc+iAfcVfc e Hk+N 

and ip^ is an eigenvector of the operator A£ +Ar A k+ ^ with the eigenvalue equal to 

k+N-l 

(5.39) = V- 

i=k-i 

For example, if we put 

(5-40) M*) = b j^M " ir^< 

Boxe q x i«oe 9 

then 

(5.41) Mx) = (1 _ g)2 ^ • 

Examples given above show the importance of the factorization method in the 
theory of the orthogonal polynomials and quantum mechanical spectral problems. 
In [5J the wide class of g-Schrodinger equations, which are solved by this method, 
is presented. 

5.4. The case gk = const and r is a fractional map. This example is similar 
to previous one but t is kept arbitrary. We still consider the case when hk = 1 and 
gk = const, fceNU {0}. From equation 113. 9|) one has that that case is possible for 
example if 

(5.42) B k {x) = b k (x - t{x)){t-\x) - x) 
and 

(5.43) {^ k (x)) 2 r] k (x) = a k 

for some constants a k , b k . These conditions are compatible with transformation 
rules H2.11j)J2,l()j) an d CHHt by following relations 

(5.44) a k+1 = ^a k 

(5.45) b k+ i = g k b k . 

The equations l|3.9l) reduce to quadratic equations with constant coefficient and 
constants g k are given as its roots. We will treat function (po as arbitrary as well 
as constants do, bo, c k , d k . In this case the equations ll.'i.2ll assume the form 
(5.46) 



-ak 



i>k(T~(x))+(b k +a k ~\ k )ip k (x)-b k (p k (T 1 (x))(t 1 (x)-x)i/) k (T 



<Pk{x){x - t(x)) 
Let us consider the equations 
(5.47) A k ip k = 0. 

It takes the form 

(5-48) Mx) = 7 rrMr(x)) 



(X - T{x))ipk{x) 



and 



(5.49) (Mx)) 2 Pk(x) = ^d T T(x) (Mr(x))) 2 p k (T(x)). 

a-k 
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If we postulate the following form of the solution to H5.49J1 

(5.50) (Mx)) 2 pk(x) = (x- t(x))V(x), 
where A 6 Z, then we get 

(5.51) a k (x) = ^(d T T(x)) x+1 a k (T(x)). 
Necessary condition for this equation to have non-zero solution is 

(5.52) (^.) X+1 =Q TT{T oo {x) y 

Let us now consider as an example the fractional map 

/ \ ax 

(5.53) t(x) 



(a - l)x + 1 



for a > 0, a ^ 1, which preserves the interval (0, 1) and has two fixed points — 
and 1. Let us set S = O(^). One can find formula for r k : 



a k x 



(5 - 54) + 
It is easily computed that 

(5.55) (9 t t)(x) = 



(-l + a 2 )x + l 
and 

(5.56) (d T T)(T k (x)) = a- {ak ~ llr ' 



'{a k + 2 -l)o;+l 
Iteration of 115, 5 1|) gives us that 



(5.57) a k (x) 



a k (0) 0<a<l 

f <a-l)x+l\ A+1 /-.x . -, 

( / M 2 j a*(l) a>l 



For A ^ -1 the function {ip k (x)) 2 p k (x){x — t(x)) has the pole at or 1, while for 
A = —1 it is constant. Thus V'fe ^ 7~Lk for any k £ N U {0}. It is also possible to 
show that kernel of A^ +1 also is empty. 

We can still consider solutions to Q5.47H as functions and use ij3,5ll to construct 
new solutions to H5.46J1 . For example, let us fix 

(5.58) po(aO = 

X — T(Xj 

From 113. 8|l and assumptions of this example it follows that 

(5 - 59) MX)= D k G k T k {x)-T k ^xy 

where G k — g k ■ ■ ■ go and D k = d k ■ ■ ■ do- Equations II5.47|I assume the form 

k ( X ) ~T k+1 (x) 



(5.60) M*)=D k G k T 



x — t(x) 
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Its solutions can be found by iteration and are given by 

(5.6!) w ,,_gfi£_J3£^L_^ fc(p) 

j=0 ^ - 1 

for < a < 1 and (GkDk)^a = 1 or 

(5.62) ^(z) = _[_[ ^J+TZ2 M 1 ) 

for a > 1 and (GkDk)i = a. 

From these solutions we can construct new by l|3.5ll . 

Appendix A. t-Riccati equation and t-Darboux transform 

We will present here additional information about r-difference equations and 
methods of solving. Most of this appendix is a direct generalization of 
Let us begin this section by solving the equation 

(A.l) d T ip = ip 

with initial condition ip(T°°(x)) = 1, under he assumption that r has the contracting 
property \t(x) — r(y)\ ^ M\x — y\, where < M < 1. Then, by the iteration one 
finds 

(A. 2) ip(x) = exp T (x) := TT — 

y ' ^ w tW 11 1-t"(.t )+T n+1 (x) 

We shall call exp T the r-exponential function. 

Let us formulate one more identity which was used throughout the paper. 



(a.3) n ( T " i? ( a; )) = cx p if x 



n=0 

It follows from the observation that 



(x)t-T(t) 



(A.4) lnf[ ( T n F(x)) =J2T n \nF(x) 



and from the definition of r-integral ll.fljl . Let us note that exp and In appearing 
in (|A.3j) are standard exponent and logarithm. 
From ijA.aj) it follows that the solution to 

(A.5) Sri? = fi> 

is given by 

(A.6) «.) = exp (£ u ^±|M1m) *(,-(.». 
Now, we are going to investigate the functional equation of the form 

d T tp(x)\ _ r , > f tp(x)\ _ (a(x) b(x)\ I V>(x) 



(A ' 7) \dMx)) - A{X) \<p(x)J - \c(x) d(x)J \ si. )) • 

By the substitution 

(A.8) A(x) = I - (x - t(x))A(x) 
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it can be equivalently written as 

(T1>(x)\ = ( , U(x)\ fa(x) b(x)\ Mx) 



(A - 9) {t<p{x)J - MX) \<p(x)J ~ \c(x) d(r) j \ ,:ir) 

If we introduce the function 

<p(x) 



(A.io) "M-^y 

the equation l|A.9J) assumes the form 

(A ID u(t(x)) - d(x)u(a:) ± c(x) 

(A ' llj U[T[X)) ~ b(x)u(x)+a( X y 

which is equivalent to 

(A. 12) d T u(x) — c(x) + d(x)u(x) — a(x)u(r(x)) — 6(x)u(x)u(t(x)). 

Equation (|A. 121 is generalization of g-Riccati equation and thus we will call both 
ilA.llj) and llA.12j) the r-Riccati equation. 

Equation 13. If) considered in Section can be obtained from ljA.91) by putting 

/ A-/3(g) 7 (g) \ 
(A. 13) A(x) = °W J . 

The formal solution to llA,9(l is given in terms of the infinite matrix product 
(A. 14) Aoo(a;) := lim A(r n (x)) • • • A(t 2 (x))A(t(.t))A(x) 

n — >oo 

by the formula 

(A 15) W-AW- 1 ^ 00 ^ 



Sometimes we are able to compute resolvent function l|A.14j) explicitly. 
Proposition A.l. 

If a, b and d are continuous at t°°(x) and c(x) = then 

(A,e, u,).^* F <*> \ 



V ^P/r-MF^W/ 



where F(x) is given by 
(A.17) 



^, ..„(/■ ^u*)f ««■,[/' ^fU K < 

W r oc (;r) < - r(f) / J T ^ {x) (t - T{t))a(t) \ J T ^ [x) s - t(s) 



Proof: From the fact that A is upper-triangular we conclude that Aoo is given 

by 

^ *-M-(-o W 

where a^x) = n^Lo 1 "" ^) and rf °o(x) = Yl^ =0 T n d(x). Since 

(A.19) A 00 (x) = A DO (r(x))A(x) 

we see that F(x) satisfies the equation 

(A.20) F(x) = rf(x)F(r(x)) + &(x)d 00 (r(x)). 
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Solving IIA.2(1jl we get 

From (IA.3I . 1IA.6I) and <Q1) we obtain now l|A.17jl. □ 

We should note that if one considers equations ijA.llj) or (j A . 1 2ll then the matrix- 
valued function A is defined up to a function factor and this observation can be 
used to regularize the product l!A.14j) . 

We are going to consider the following gauge-like transformation 

(A. 22) A(x) — ► A'(a;) = D(t(x))^ 1 A(x)D(x), 

for some D : X — » GL(2, C), which transforms the functional equation (lA.Qji into 
the one of the same type with matrix A'. Its solutions are given by 

(A.23) (?Sfl^Dto-iM4 



The resolvent of new equation is given by 

(A.24) AooOk) — > D(0)- 1 A oo (x)£>(x). 

We call this transform the r-Darboux transform. It is a generalization of q- 
Darboux transform introduced in ^3]. Sometimes we are able to reduce A to such 
A' that we can compute l|A.14|) explicitly. This is the case, for example, if A' is 
upper-triangular. Applying this method, we can obtain the result for r-Riccati 
case, which generalize this of g-Riccati equation of [13] . 

This is a classical result for the Ricatti equation generalized to q-case in JS| and 
here it is presented in r- version. 

Proposition A. 2. 

Let uq be a particular solution of l|A.ll|l . The general solution of (fA.Qf) is then 
given by the formula 

= ma( t ) + mMt) d \ f r b {t ) 




t-r(t) J y J TOO{x) a(t) + b(t)u (t) 

a(s)+b(s)uo(s) 

(A.25) <- M , [ / d T s } d T l + A 



(A.26) V {x)=u {x)^{x)+Be^\-I ^-WM^)) +d(t) ^ 




where the constants A and B are related to the initial conditions in the following 
way 

(A.27) A = iP(t°°(x)) 

and 

(A.28) B = -1>(t°°(x))uo(t 00 (x)) + v(t°°(x)). 
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Solution of l|A.llj) is then given by 
(A.29) 

u\x) = u (x) + - 



/or any f £ 1. 

Proof: In order to obtain the general solution of l|A.7ll . we apply the r-Darboux 
transform l|A.22jl with D{x) of the form 

1 X 

Since u$ is a particular solution of IIA.llj) . the transformed matrix A' (a;) is upper- 
triangular. Solution given by IIA.2J) and l|A.26jl is obtained by substituting into 
i|A.lfiJI the formulae IIA.lfill and 1IA.17JI from Proposition IA.1I and transforming 
formula l|A.24j) for resolvent. Formula l|A.29j) follows from IIA.lOjl . l!A.2jl and l|A.26jl 



(A.30) D(x) 



by putting t := -r. □ 
We define the transform Bf by 

(A.31) Btu := u l . 

It is easy to check that Bf , t S R, form a one-parameter group 

(A.32) BtoB+ = B+ +t , 

Moreover the solutions u' 1 (a:) , u* 2 (x) , u ta (x) and u' 4 (a;) do satisfy the unharmonical 
superposition principle: 

33 (u^(x)-u t Hx))(u^(x)-u t -(x)) = (U-t 3 )(h-t 2 ) 

The proof of (IA.32I and (IA.33I) is straightforward. 

Let us mention another application of r-Darboux transform. When A(r°°(a;)) = 
I we call IIA.QJl the regular equation. Solutions of regular equation are finite and 
nonzero at t°°(x). Let us introduce singularity (or zero) at this point by the r- 
Darboux transform (IA.22I) - (I A. 231) with 



'{X -T°°{x)) Sl 

(X-T°°(x)) 

Thus one has 



(A.34) D(x) — i , ^oc^5 2 



( - a'i .c i l^w-^wJ 

C W ( r(3; )_ r oo (;r)) o 2 a^j ^ T(a; )_ T oo (2 . ) J 



and the solution to (|A. 1 If) with A'(x) is 

(A.36) u\x) = {x- t°°(x)) Si - S2 u{x), 

where u(x) is a solution to l|A.llj) with A(x). 

Hence one can have u'(t°°(x)) — oo for 5\ — 5 2 < while u(t°°(x)) < oo. Thus 
we can apply the transformation l!A.34|l to regularize the singular equation of the 
type (IA.35II to the regular one. It was applied implicitly for example in H5.50I) . 
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Infinite products in previous formulae were always assumed to be convergent. 
We will give now a condition for this to be the case. 

Proposition A. 3. 

Let the following conditions be satisfied: 

i) jy^Lo \ rk ( x ) ~ T k+1 (x)\ < oo, e.g. it has a place if t is a contracting map; 

ii) A is continuous atT°°(x). 

Then product (I A . 14|) is convergent. 

Proof: We will show that under these assumptions P n = A(r n (x)) ■ ■ ■ A(x) is a 
Cauchy sequence and thus convergent. We can estimate the norm of the difference 
between two terms of this series by 



(A.37) 



Pn - Pn+k\\ < A(r n+ k (x)) • ■ • A(r" +i (x)) - 1 ||P„| 



Let us show that it can be made arbitrarily small. To that end we show first that 
\\P n \\ is bounded. 



< oo. 



(A.38) ||P„|| ^ H(l + \r k {x) - r k+1 (x)\ ||A(r fc (x))||). 

fc=0 

This product is convergent if and only if 

oo 

(A.39) ^2\T n {x)-T n+1 {x)\ A(r n (x)) 

71 = 

Assumptions of the proposition guarantees that this is the case. Hence 

oo 

(A.40) ||P n || ^ + \t\x) - T k+ \x)\ ||A(r fc (x))||) < oo. 

k=0 

Now we show that for n big enough || A(r™ +fe (a;)) • • • A(t" +1 (.t)) — l|| is arbitrarily 
small. We have the inequality 



(A.41) 



A(T n+k (x)) ■ ■ ■ A(T n+L (x)) - 111 + 1 < 



S$ (1+ \T n+ *{x) -T n+k+1 (x)\ A{T n+k (x)) )---(l+\T n+l (x)-T n+2 (x)\ A(r" +i (.T)) 



n+k 



Since IIA.39J1 hold then l!A.41|l tends to 1 as n — > oo. Thus we have completed the 
proof that P n is Cauchy sequence. □ 
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